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Abstract. We analyze the dynamics of various kinds of correlations present between 
two initially entangled independent qubits, each one subject to a local phase noisy 
laser. We give explicit expressions of the relevant quantifiers of correlations for the 
general case of single-qubit unital evolution, which includes the case of a phase noisy 
laser. Although the light field is treated as classical, we find that this model can 
describe revivals of quantum correlations. Two different dynamical regimes of decay 
of correlations occur, a Markovian one (exponential decay) and a non-Markovian one 
(oscillatory decay with revivals) depending on the values of system parameters. In 
particular, in the non-Markovian regime, quantum correlations quantified by quantum 
discord show an oscillatory decay faster than that of classical correlations. Moreover, 
there are time regions where nonzero discord is present while entanglement is zero. 
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1. Introduction 

Total correlations present in a system can be separated in a purely quantum part 
and a classical part. Quantum correlations can be classified as those associated with 
non-separability (entanglement) and other quantum correlations, together quantified 
by the quantum discord (QD) [H [2]. For quantum information processing it is 
relevant to know the dynamics of these different kinds of correlations. Entanglement 
dynamics can be considered well understood, in its general lines, for bipartite quantum 
systems interacting with quantum environments (independent or common), presenting 
phenomena like sudden death [3], revivals [H [5l E] or trapping [7], depending on the 
Markovian or non-Markovian nature of the environments. Dynamics of QD has also 
been investigated for two-qubit systems in the presence of both Markovian [HI |9] and 
non-Markovian [101 [H] quantum environments. 

In this paper we analyze the dynamics of correlations in a two-qubit system where 
each qubit is subject to a phase noisy laser modeled as a classical field. This is a 
commonly used model for the interaction between light and matter, and thus the 
dynamics of correlations for this situation merits investigation. The broader question 
underlying our work is whether all kinds of behavior of correlations that can be described 
using a quantum environment also can be well described using a classical approximation 
of the environment. One might for example expect that a classical environment should 
not be able to store quantum correlations on its own, and that therefore revivals of 
quantum correlations in the system may be affected. For the phase-noisy laser, are 
there qualitative features that will be lost by treating the field classically as opposed to 
quantum-mechanically? Among other things, we find that this model can describe both 
decay and oscillatory behavior of quantum correlations. 

2. Model: a qubit subject to a phase-noisy laser 

Our system consists of a pair of qubits (two-level atoms), A and B, each driven by a 
local phase noisy laser. Each atom under the action of its own laser is described, in a 
rotating frame and for resonant atom-field interaction, by the Hamiltonian [12] 

H = X [(7„e^*W + o-+e-^*W] , (1) 

where the laser is described as a classical field with a randomly fiuctuating phase $(t). 
The interaction between each qubit and its local field mode is assumed to be strong 
enough so that, for sufficiently long times, the dissipation effects of the vacuum radiation 
modes on the qubit dynamics can be neglected. This could be feasible by considering, 
as a qubit, an atom in a cavity subject to a resonant interaction with the phase noisy 
laser but out of resonance with cavity mode frequencies in order to inhibit effects like 
spontaneous emission. In this phase-noisy model the phase undergoes a Wiener process, 
i.e. $(t) is white noise with a correlation function ($(t)$(t -|- r)) = 2c/5(r) where 
(i is a diffusion rate. In Eq. ([1]), A is the atom-field coupling strength, a+ = |1)(0| 
and cr_ = |0)(1| are the atomic raising and lowering operators, where |0) and |1) are 
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the ground and excited state of the atom, respectively. The field correlation function, 
corresponding to the above phase correlation function, is a complex colored noise e***-*^ 
described by (e^*(t)e-^*(t + r)) = e"'^^. 

The Hamiltonian of Eq. ([1]) leads to a local-in-time non-Markovian master equation 
of the form [12] 

3 

p{t) = '^ili(^ip(^i - lip), (2) 

i=l 

where (jj (z = 1, 2, 3) are the atomic pseudo-spin Pauli matrices and 

This master equation is derived by the Nakajima-Zwanzig projection operator method 
and the functions ri,r2 can be expressed in terms of the system parameters. In 
particular ri(t) = e'f [cosh {\ts/d? - 16A2) + c/sinh {\t^Jd? - IGA^) /yjd? - IGA^J , 
while the expression for r2(t) is cumbersome and is not given here. Because the Fj 
are time-dependent, the master equation of Eqs. ([2]) and ([3]) has a quasi-Lindblad form, 
that is a master equation that resembles the Lindblad form but has time-dependent 
decay rates 7j (which also can be negative) [12]. Moreover, the above single-qubit 
master equation is a unital master equation [13]. A master equation dp/dt = Ct{p) is 
defined to be unital if the maximally mixed state |J is a fixed point, that is if jCt{I) = 0. 
Unital master equations can be solved by Kraus-type decomposition methods giving the 
single-atom reduced density matrix elements at time t as 



pii(t) = [(1 + A3)pii(0) + (1 - A3)poo(0)]/2, 

pio(t) = [(Ai + A2)pio(0) + (Ai - A2)poi(0)]/2, (4) 

with pii = 1 — Poo, pio = Poi- Aj(t) = Ai are time-dependent and related to the 

rates 7j by 

A,(t) =e-2/o<i*'[^.(*')+^^(*')l, (5) 

with the conditions A^ -|- Aj < 1 + A^, where {i,j, k} run over the cyclic permutations 
of {1,2,3}. Note also that Aj{0) = 1 and Aj{t) > 0. Using Eq. in Eq. (ED we 
have Ai = A2 = r2, A3 = Fi and the resulting reduced density matrix elements 
for the single atom driven by a phase noisy laser reduces, from Eq. (jlj), to pu{t) = 
|[(1 -|- Fi)pii(O) -|- (1 — Fi)poo(O)] and pio = 2F2Pio(0). It has been shown that for 
d/X < 4, Fi oscillates while F2 shows oscillatory behavior only for d/X < 2.606 [12]; for 
d/X > 4, neither Fi nor F2 shows oscillatory behavior. Finally, in the limit d/X ^ 0, the 
two decay rates Fi, F2 become periodic functions of At, more specifically Fi — )• cos(2At) 
and F2 — 7- cos^(At). In this limit, the probability distribution of the phase is a steady 
state distribution, so that the effective dynamics is the statistical average, equally 
weighted, of all the evolutions with phase between and 27r. An analogous but simpler 
setup, described by a uniform phase distribution with only two values (0,7r), has been 
investigated in Ref. [E 
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3. Quantifiers of two-qubit correlations 

To obtain the expressions of the correlation quantifiers we need the two-qubit density 
matrix elements. We construct the two-qubit reduced density matrix at time t by the 
knowledge of the evolution of single-qubit reduced density matrices, according to a 
standard procedure [1], with the single-qubit reduced density matrix evolution given 
by Eqs. (jlj) and (Q. Thus, we obtain the explicit expressions of the two-qubit density 
matrix elements at time t in terms of the functions Aj for any initial two-qubit state. 
These expressions are however quite cumbersome and are not reported here. 

We take as two-qubit initial states the extended Werner-like (EWL) states [5] 

p* = r|<|.) ($1 + 1^/4, f = r\^){^\ + ^h, (6) 

where r indicates the purity of the initial states, I4 is the 4x4 identity matrix, 
1$) = ci|01) + (3e'^\10) and |^) = a|00) + I3e'^\ll) are the Bell-like states where are 
non- negative real numbers and a"^ + (5"^ = 1. These are mixed states reducing to Werner 
states for a = (3 = ±1/^/2 or to Bell-like states for r = 1. The density matrix elements 
of the EWL states are such that the resulting density matrix has an "X" structure with 
nonzero elements only along the main diagonal and anti-diagonal. 

In the standard basis B = {|1) = |11), |2) = |10), |3) = |01), |4) = |00)} and for 
general different environments characterized by different values of Af (z = 1,2,3; S = 
A,B), the time-dependent two-qubit density matrix elements for /)* are 

PJW = ^ {1 - ^ [AX + (-1)^(1 - 2«')(A^ - Af )] } , 

Pz!W = \{l+r [aX + (-l)'(l - 2a^)(A^ + Af )] } , 

pt,,_,it) = ^ [/(A) cos5 + i{-l)'^'9{A) sin 5] , (7) 

where j = 1,4, I = 2,3, k = 0, 1, /(A) = A^ + A^ and ^(A) = A^ + A^- The 
density matrix elements for the initial state p* are obtained from Eq. ([7]) by changing 
1^2 and 3 -H- 4. Note that under the dynamical conditions here considered, the 
two-qubit state maintains an X structure during the evolution. 

In order to describe the entanglement dynamics we use the concurrence, which 
for an X state is given by ^ Cf{t) = 2max{0, Ki{t), K2{t)} where Ki{t) = 
\p23it)\- ^Jpii{t)pAi{t) and K2{t) = |pi4(t)| - VP22(t)p33(t). The EWL states of Eq. (El) 
present the same initial value of the concurrences, C(0) = 2max{0, (a/3 -|- l/4)r — 1/4}, 
from which one finds that there is initial entanglement when r > r* = {1 + 4a/3)~^. 
Using the time-dependent density matrix elements of Eq. ([7]), it is readily seen that 
in our system the concurrence at time t is the same for both the EWL states 
of Eq. ([6]), that is C*(t) = Cp{t) = C(t). For example, for initial Bell states 
{r = l,a = 1/a/2, 5 = 0, vr) and different local conditions, concurrence is given by 
Cit) = I {A^Af + A^Af + A^Af - 1). 

In order to quantify total correlations, T, present in the two-qubit system and 
to distinguish a quantum V and a classical part J7 of them, we use the notion of 
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quantum discord [UlZj- The calculation of discord and classical correlations requires a 
maximization procedure and this has been analytically solved only for certain class 
of quantum states (Bell-diagonal and X states) [Ml [IT]. In particular, here we 
give the explicit expressions of these quantifiers for the initial condition a = (3 = 
l/\/2,S = 0,7r for which the two-qubit density matrix always has a Bell-diagonal form 
Pb = [I (8) I + X]j=i CjCTi ® o"j]/4. For this class of states, the following expressions for T 
and J, with V = T - J, hold [T6] : 

r=2 + Y,^^og\l J = Y.^-^^^\og[l + {-!)% (8) 

i,s i 

where z = 1, 2; s = ±, c = max{|ci|, |c2|, jcsl}, = (1 ± Ci ± C2 — C3)/4, = (1 ± Ci =F 
C2 + C3)/4. If the initial state is /)*, the q coefficients are ci = rmax{A^Af , A^A^}, 
C2 = rmin{A^Af , A^A^} and C3 = — rA^A^. If the initial state is p*, previous 
coefficients change as ci — )■ ci, C2 — > — C2 and C3 — — C3, that is a relabeling of the 
A eigenvalues: the quantifiers T, V and JT" thus coincide for both initial states. 

4. Dynamics of correlations due to a phase-noisy laser 

We use the general results of the previous section to analyze the dynamics of correlations. 
In particular, for initial EWL states in identical environments {Vf = V^ = Fj, j = 1, 2) 
we obtain C*(t) = C^(t) = C = \ {Aral3V\ + rV\ — 1). In Fig.[T]^a) we plot concurrence 
evolution for initial Bell states, for two different values of the ratio d/X = 5,0.1. The 
behavior for more general initial mixed states can be shown to be qualitatively similar. 
One sees that while concurrence presents Markovian-like decay for d/A = 5, it is subject 
to non-Markovian revival for d/X = 0.1, as already found in other systems [H [T3]. In 
order to evidence the fact that for large times the two-qubit system goes toward a 
maximally mixed state as a consequence of the unital evolution of the single qubits, in 
Fig. Wljo) we plot the evolution of the von Neumann entropy, S = — Tr{p(t) logp(t)}, for 
the same values of d/X (von Neumann entropy is the same for both initial EWL states 
^pi't) = Sp(t) = S). It is seen that, for any value of d/X, S tends to its maximum value 
2, corresponding to the two-qubit maximally mixed state, with oscillating behavior when 
memory effects are present {d/X = 0.1). 

We can also separately investigate the evolution of quantum and classical 
correlations. In Fig. |2]we plot the evolution of total, quantum and classical correlations 
of Eq. (|8]) for the value of the ratio d/X = 0.1 for initial Bell states. From the plot one 
sees that, in this non-Markovian case, total and classical correlations show an oscillating 
decay while quantum discord has a non-purely oscillating decay (i.e., there are points 
of discontinuity for the derivative during the time regions where discord is very close to 
zero) and it decays much faster than classical correlations. On the other hand, as seen 
from the inset of Fig. |2l discord definitively disappears later than entanglement, with 
time regions when it is different from zero while entanglement is absent. 
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Figure 1. Time evolution of concurrence C (panel (a)) and von Neumann entropy S 
(panel (b)) starting from initial Bell states {r = l,a — l/\/2,S = 0, tt) for the ratios 
d/X = 5 (black solid curve) and d/\ = 0.1 (blue dashed curve). 
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Figure 2. Time evolution of quantum discord V (black solid curve), classical 
correlations J (red dashed curve) and total correlations T (blue dot-dashed curve), 
starting from initial Bell states {r — l,a ~ l/^/2^5 — 0,7r) for the d/X ~ 0.1. In the 
inset V (black solid curve) and concurrence C (blue dashed curve) are plotted. 
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5. Conclusions 

In this paper we have analyzed the dynamics of correlations between two initially 
entangled independent qubits each locally subject to a phase noisy laser. We found 
explicit expressions for various quantifiers of correlations, such as concurrence and 
quantum discord, in terms of the decay rates present in the single-qubit unital master 
equation which here describes the action of phase noisy laser. We discussed how the 
dynamics of correlations depends on the ratio between the phase diffusion rate d and 
the atom- laser coupling strength A. 

Although the light field is treated as classical, this model can describe revivals 
of quantum correlations. For large values of dj\ Markovian-like decay occurs, while 
non-Markovian effects become relevant for small values of this ratio, giving place for 
example to revival of entanglement and of discord. Moreover, in the non-Markovian 
regime, quantum discord presents an oscillatory decay faster than that of classical 
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correlations, with time regions where it is nonzero and entanglement is zero. It is 
therefore evident that the phase diffusion rate and coupling strength of each atom-laser 
affect the qualitative time-behavior of quantum correlations in a crucial way. 
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